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Abstract 

We discuss the role of integrated chance constraints (ICC) as quantitative risk constraints in as¬ 
set and liability management (ALM) for pension funds. We define two types of ICC: the one period 
integrated chance constraint (OICC) and the multiperiod integrated chance constraint (MICC). As 
their names suggest, the OICC covers only one period whereas several periods are taken into account 
with the MICC. A multistage stochastic linear programming model is therefore developed for this 
purpose and a special mention is paid to the modeling of the MICC. 

Based on a numerical example, we firstly analyse the effects of the OICC and the MICC on the 
optimal decisions (asset allocation and contribution rate) of a pension fund. By definition, the 
MICC is more restrictive and safer compared to the OICC. Secondly, we quantify this MICC safety 
increase. The results show that although the optimal decisions from the OICC and the MICC differ, 
the total costs are very close, showing that the MICC is definitely a better approach since it is more 
prudent. 
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1 Introduction 


A pension fund is any plan, fund or scheme, established by a company, governmental institution or 
labour union, which provides retirement incomes. The actuarial present value of current and future 
payments constitutes the total liability of the fund. The pension fund receives contributions from 
its active members and/or the employer. This money (considered as the total wealth or total asset) 
is invested in a wide range of assets. The asset allocation is made in such a way that it guarantees, 
to a certain extent, the payments of future obligations. That is not so trivial: assets yield random 
returns and future benefits are not known with certainty. An asset liability management study 
(ALM) provides a rich theoretical background to address that issue. Its goal is to determine the 
adequate asset allocations and contribution rates in order to guarantee the payment of current and 
future pensions. 

The use of ALM methods has a long tradition in pension funds. At the beginning, it has started 
with deterministic methods. According to those methods, the future cash flows are estimated and 
assumed to be certain; the wealth is mainly allocated to bonds considered as risk free. Bonds are 
chosen in such a way that their related incomes correspond to yearly pension payments. These 
models are essentially based on immunization and cash-flow matching; see, for example, Koopmans 
[18] and Redington m- Deterministic methods have proven to be inefficient since uncertainty turns 
out to be increasingly difficult to handle. The problem of uncertainty is somehow taken into account 
by stochastic methods, especially by the way of surplus optimization theory. This approach is often 
based on the efficient portfolio theory of Markowitz [25] . The literature has flourished in that 
held and we can cite Sharp and Tint [35] and Leibowitz [ 22] among others. However, the pension 
fund problem is a long term problem with a horizon span of approximatively 30 years. Hence, its 
model should be dynamic. Furthermore, regulations often impose many types of constraints. Those 
matters are hardly taken into account by surplus optimization methods. In practice, simulation 
methods are commonly used due to their ability to incorporate the above issues. Initially, they 
consisted on defining a set of feasible allocations and contribution rates, and chosing the best one 
in some sense. The choice is based on the simulation of the future paths. Due to the technical 
innovations, these methods have significantly evolved with the work of Wilkie |39] and Ahlgrim and 
al. [I] concerning economic scenario generation. Mpller and Steffensen [26] provide different tools 
for valuing the pension fund liabilities. Recent years have also seen the emergence of methods known 
as stochastic programming. 

Often based on simulations due to its complexity, stochastic programming gives a flexible and 
powerful tool for ALMs. Its importance lies in its ability to bring together many kinds of features 
in a common framework. Moreover, assets and liabilities are all influenced by many sources of risk 
and the risk aversion is accommodated; the framework has a long time horizon split into subperiods 
(multistage ); the portfolio can be rebalanced dynamically at the beginning of each subperiod; all 
these are incorporated in a single and consistent structure while satisfying operational or regulatory 
restrictions and policy requirements. Multistage stochastic programs (MSP) models have been 
applied to ALM for pension funds by Carino and al. [7], Consigli and Dempster [9j, Kusy and 
Ziemba m and Kouwenberg and Zenios |19| . 

The ALM model in this paper is a MSP, for which, we minimize the total funding cost under risk, 
legal, budget, regulatory and operating constraints. The total funding cost is composed of regular 
and remedial contributions. Regular contribution constitutes a certain proportion (contribution 
rate) of the total salary whereas remedial contribution is an additional financial support provided 
by the employer (or a sponsor) whenever the solvency target is in question. More specifically, we 
focus on the risk constraints, which are of integrated chance constraints (ICC) type in this work. 
As an alternative to chance constraints (CC), ICC is computationaly of great interest; in particular 
when a quantitative risk measure is preferable. We define the funding ratio as the ratio of total 
asset over total liability. Our goal is to meet a certain funding ratio, called here target funding ratio, 
at the end of each subperiod. For a predefined target funding ratio, the ICC put an upper bound 
on the expected shortfall, i.e. the expected amount by which the goal is not attained. Haneveld 
and al. m and Drijver m pionnered the application of ICC in ALM for pension fund. However, 
the risk parameter considered in their models is neither scale free, nor time dependent. Our model 
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is close to Haneveld and al. m with the particularity that the risk parameter is a linear function 
of total liability. Then, it becomes unvariant with respect to the size of the fund as well as time 
dependent. 

We define two types of ICC: the one period integrated chance constraint (OICC) and the multiperiod 
integrated chance constraint (MICC). As their names suggest, the OICC covers only one period 
whereas several periods are taken into account with the MICC. A multistage stochastic linear 
program is therefore developed for this purpose and a special mention is paid to the modeling of 
the MICC. 

The rest of the paper is organized as follows. In section [2] the theoretical background, the dynamics 
and the ALM optimization problem are extensively detailed. Section [3] defines the risk contraints 
and shows how CC leads to ICC. Moreover, OICC and MICC are introduced and their stochastic 
linear program reformulations are derived as well. In section [dj a numerical example is examined 
from the perspective of a defined benefit fund that invests in stocks, real estate, bonds, deposits 
and cash. All numerical results are implemented using the solver CPLEX in the mathematical 
programming language AMPL. We first analyse the effect of the risk parameter on the optimal 
decisions. This section finishes by a brief comparison of the two ICC. Section [5] concludes the paper. 


2 Settings 

This chapter introduces the dynamics of the ALM model with its specific features. 

2.1 Multistage recourse models 

In this section, we describe the classical architecture of the multiperiod decision framework. The 
model’s setup presented here ressembles mostly to Haneveld and al. [fT3] , 

Since we aim for strategic decisions, we model the ALM process over a number of years and one 
set of decisions is taken each year. We discretize time accordingly so that the model has a (finite) 
number of one-year time periods. Consequently, we assume that the ALM model has a horizon of 
T years from now, split in T subperiods of one year each. The resulting years are denoted by an 
index t, where time t = 0 is the current time. By year t (t = 1, • • • , T), we mean the span of time 
[t— 1 ,t). We define 

% := {t,t + 1, • • • , T} . 

We assume that uncertain parameters (e.g. asset returns) can be modeled as random variables 
with known distributions. At each time t G To, the pension fund is allowed to make decisions 
(corresponding to yearly corrections), based on the actual knowledge of parameters. During each 
one-year period, a realization of the corresponding random parameters becomes known (e.g. assets 
return during that year). That is, the concept underlying our model is the following sequence of 
decisions and observations: 

decide observe decide observe decide observe 

X 0 -—> X | (U 7 1 1 Xn 1 ~—> lot 

where Xf is the vector of decision variables at time t G 7o, and vector uit, t G 71 models all economic 
events which are the source of uncertainty and risk for the pension fund management, which, in our 
case, are asset returns as well as random contributions and liability streams. Time t is assumed 
to be the end of the financial year t. We assume that a financial year coincides with a calendar 
year. At time t G 7o, decisions Xf are taken with full knowledge of the past [0, t] but with only 
probabilistic informations about the future (t, T]. 

Uncertainty in the model is expressed through a finite number S of sample paths spanning from 
t = 0 until t = T called scenarios. That is, we assume the random variable follows a discrete 
distribution with S possible outcomes. Each scenario represents a sequence of possible realizations 
of all uncertain parameters in the model. As explained above, uif is the stochastic vector process 
whose values are revealed in year t. Then, the set of all scenarios is the set of all realizations 
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t — 0 event t — 1 event £ = 2 event t = 3 
Figure 2.1.1: A scenario tree with 40 scenarios and 66 nodes. 


uj s := (cjf,--- ,w|),s £ 5 := {1, • • • , 5} of u> := (uq,--- ,ur). Scenario s has a probability p s , 
where p s > 0 and ^ 5=1 p s = 1. It represents a description of possible future, starting just after 
t = 0. If we assume that we can observe the "state of the world" at time t, (0 < t < T ), then 
there is a unique history of realizations of (wi, • • • ,0Jt- 1 ) leading to that state, but the future as 
seen from time t may unfold in several ways. That is, there are several distinct scenarios which 
share a common history up to time t. A suitable representation of the set of scenarios is given by 


a scenario tree (see Figure 2.1.1). In respect to Figure 2.1.1 we define the node as the possible 


outcome of the stochastic vector ut at a given time t € To- Each path of ut from t = 0 to t = 3 
represents one scenario; each node of the scenario tree has multiple sucessors, in order to model the 
process of information being revealed progressively through time. By convention, the scenarios are 
numbered top-down by their end node. The arcs in the tree denote realizations in one time period. 
We assume here that, for a specific decision time t E To, the numbers of realizations in one time 
period descending from the current nodes are identical. 


For example in Figure 2.1.1 we have a 3-year horizon scenario tree with 40 scenarios. Over the 
first period starting from time 0 to time 1, there are five possible realizations. From each of these 
realisations, we have four possible outcomes over the second year; each of them is a conditionnal 
realisation as it depends on the preceding node. Over the third period, each of the second period 
observations can lead to two possible outcomes. All this gives a branching structure of 1 — 5 — 4 — 2 
and leading to a total ofS' = 5x4x2 = 40 possible scenarios. 

A multistage recourse model is an optimization problem defined on such a scenario tree. Consid¬ 
ering the remaining future represented by the subtree rooted at (t,s), optimal decisions are taken 
for each node (t,s) of the event tree, given the informations available at that point. Optimality is 
defined in terms of current costs plus expected future costs, which are computed with respect to the 
appropriate conditional distributions, Vlerk and al. m 

Ideally, one would like to make different decisions for every path at every t E To, but this would 
lead to undesirable anticipativity in the model. The simplest way to avoid this is to make one single 
decision at each time t for all paths by adding explicit constraints. That is, for any two different 
scenarios si and s 2 (si, S 2 E S and si ^ s 2 ) having the same history up to time t E To, we enforce 
A t si = A/ 2 , where Xf is the decision Xt under scenario s. For example, at the empty circle of 


Figure 2.1.1 X{ = Xf = ■ ■ ■ = Xf. 
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2.2 Dynamics 

2.2.1 Assets 

In this paper, we are considering a buy and hold model applied to a DB plan in which one seeks 
to minimize the expected cost of funding. In this respect, dynamics for both assets and liabilities 
should clearly be specified. 

At initial time t = 0, the exact levels of wealth and liability are available to the decision maker who 
has to decide, each period, how to rearrange his portfolio in order to cover liabilities and, at the 
same time, to achieve high returns on the financial market. The higher the returns are, the lower 
the contribution rate could be. Let denote At the total amount of wealth at time t E To- The total 
wealth is allocated into d classes of assets and in cash. Let k € 1C := {1,..., d} denote the asset class 
index. At each decision time t E 7o, a specified amount of H^t is allocated to asset k and Ct is the 
cash amount. We can write 

d 

At = Hk,t + Ct- 

k =1 

Through buying and selling, the investor restructures his portfolio at each time t. Once the t th stage 
decision is made, the holdings H^ t can be calculated. The shares in the portfolio are then kept 
constant till the next decision time. The value of H^t is affected by the returns on the market. Let 
define C,k,t, := 1 + rk,t where rk,t is the random rate of return on asset class k over year t. 

Over year t. the pension fund pays benefits to its non-active members and receives contributions 
from its active members or/and the employer (also called the sponsor). Benefits regroup pensions 
which are paid to retirees, disability and death annuities or lump sum, whereas contributions are 
composed of yearly payments from all the active members and/or sponsor to the plan. When it 
appears that the plan is unfunded according to its solvency target, the sponsor may finance the 
deficit. As in Vlerk and al. IS], we name this funding here as remedial contribution. We then 
assume that whenever the solvency target is not fulfilled, a remedial contribution in cash can be 
obtained from the sponsor. In practice, it does not really work that way. For example in Vlerk and 
al. [.33, the remedial is only provided after two consecutive periods of underfunding. We will see 
in the model description that the parameters are set such that the remedial contribution variable 
is non-zero only under some conditions. In general for DB plans, future benefits and liabilities 
depend on company policy regulation and can be estimated whereas yearly contribution is defined 
as a certain proportion of the yearly salary. Asset allocation and contribution rate are defined with 
respect to the level of future benefits and liabilities (e.g. Switzerland). Kim [16] provides a rich 
source of informations concerning different types of pension plans and features. During year t, let 
Bent and Wt denote, respectively, the total amount of benefits paid and the level of salary. The 
variable cry is the decided contribution rate for year t + 1. For returns and cash-flow variables, 
index t means that payments occur over year t but cash-flows are accounted at the end of year. 
Accordingly, the total asset dynamic is modeled as 

d d 

At = Hk,t—i£k,t + C t —i (1 + rf) + cr't -1 Wt — Bent + %t = Hk,t + Ct, (2.2.1) 

k= 1 k =1 

for t Eli, where r/ is the risk free interest rate and Zt is the remedial contribution at time t. Before 
receiving the remedial contribution at time t, the total wealth is defined as A* and one can write 

d 

■Af = Hk,t.-i£k,t + Ct —i (1 + rf) + crt-iWt — Bent = At — Zt- (2.2.2) 

k=i 

Asset returns {Ct)t=\ := (£i,* ■■■ €k,t ■■■ id,t) t _ v pension payments Bent and salary Wt are 

modeled as stochastic processes on a filtered probability space T, (J r )^ =1 , . Obviously, at each 

decision time, At is a random variable whose distribution depends, on a first hand, on £t, kLt and 
Bent, and on the second hand, on asset allocations before t. At a specific date t, the variable At is 
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known as it can be observed. The initial wealth is defined by Aq and known at initial time t = 0. 


According to (2.2.1), total wealth At -1 at time i — 1 is allocated into the d classes of assets and 


cash. Each asset k generates a return £k,t over period [t — 1, t]. The initial wealth plus accumulated 
interest at the end of period will be augmented by the balance of external flows: contributions minus 
pension payments. This latter can be either positive or negative depending on the difference between 
contributions and benefits paid. A negative balance could be due to the fact that a company is no 
more hiring new employees. This may happen for various reasons: runoff, economic difficulties, etc. 
Obviously, the total contributions crt-iWt will decrease considerably as the total salaries decrease 
whereas Ben^ will tend to increase as people leave the fund. When at time t the total asset can not 
fulfill the pension fund solvency target, it may obtain a remedial contribution Zf. 

In order to be as close as possible to realities on the financial market, one has to consider costs of trad¬ 


ing activities. Therefore, we include proportional transaction costs c B := [c\ 


tiB 


7lB 


and c s : = (cf ... c B 
costs will lead to some changes in asset dynamics. Thus, (|2. 2. lh is then replaced by 


) 


crj) for purchases and sales, respectively. Inclusion of transaction 


At = ^ #T-i,fc£r,fc + Ct-i (1 + rf) + ctt-iWt — Ben-r = Aj 


(2.2.3) 


k =1 


over period [T — 1, T] and when t E T\ \ {T}, 


At — 'y ] Hk,t-i£k,t + C t -i (1 + rf) + crt-iWt — Bent + Zt — ( c B Bkj + c B Sk t t) 


k =1 


k =1 


— £,tHt—l + Ct—i (1 + t f) + crt-iWt — Bent + Zt — (c B Bt + P St) 
= A* + Z t - ( c B B t + c s S t ) 

= e • Ht + Ct 


(2.2.4) 


where e := (l 1 • • • l) is a (1 x d) vector. Vectors Ht := (-Hgt 


H, 


k,t 


{Bi >t 


Bk,t 


B d ,t) T and S t := {S 


H d , t ) , B t : = 


l ,t 


S k 


Sd,t) of dimension (d X 1) each, 


are, respectively, amount of asset hold, bought and sold at each decision time t E 7o- In fact, (2.2.4) 


is obtained by substracting transaction costs in the first equality of (2.2.1). At time T, no more asset 


is bought or sold: Bt = St = 0; the value of the portfolio is determined by adding all values of assets 
including the last period returns and external flows. This justifies why there is no transaction cost 


in (2.2.3). The reader should notice that variables crt , Zt, Bt, St and Ht are all decision variables. 
We denote by H^, the initial holding in asset k, k E K, and H := ( H\ ■ ■ ■ Hf ■ ■ ■ H d ) T is a 
d x 1 vector. C$ is the initial cash amount. The first stage asset allocation is determined by 


H 0 = H + B 0 - So 


with total asset 

A 0 = e • H + Co + Zq — (c b B 0 + c B So) = A-o + Zq — {c B B q + c B So) = e • Hq + Co- 
For t > 1, 

H t = ZtHt-x + B t -S t 

defines the dynamic of holding assets between two consecutive decision times. For any given ( k,t ), 
whenever S^t > 0, B^ t = 0 an d vice-versa. Transaction costs also influence the cash dynamics. 
Buying an amount Xk of asset k requires Xk (l + c B ) of cash and selling the same amount of asset 
k results in Xk (l — c B ) of cash. Initially, 

Co = Co + Zo — (e + c B ) Bo + (e — cP) So 

and for t > 1, 

C t = Ct— i (1 + rf) + crt-iWt — Ben^ + Zt — (e + c B ) Bt + (e — c B ) St 
where we assume that crt-iWt, Ben^ and Zt come in cash. 
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2.2.2 Liability and external flows 

As we consider a DB plan, total liabilities is the discounted expected value of future pre-defined 
payments. At a given time t. it represents the amount the fund has to reimburse if it has to close 
at that time. Its value has to be estimated with appropriate rules taking into account actuarial 
risks, pension fund provisions, and other relevant factors for the employer’s line of business. Let Lt 
denote the total amount of liabilities at time t. 

All quantitative models considered in this paper will be applied to the planning problem of a large 
and stable pension fund. We can then assume that the fund keeps the same structure and number 
of members over the study period. Liability, contributions and benefits are therefore invariant with 
respect to actuarial risk over the period under study. Actuarial risks regroup the random events 
that affect the number of members into the fund. However, those variables are yearly indexed with 
the general increase of wages Wt- For t G 7i, we have: 

L t = L t - i (1 + w t ); W t = W t - 1 (1 + w t ) and Ben t = Ben t _i (1 + KW t ) (2.2.5) 

and their initial values Lq, Wo and Beno known at t = 0; k > 0 is a model parameter. In practice, 
the pension payments Ben^ are often indexed with a certain rate which is a function of the inflation 
rate. In order to reduce the complexity of our model, we assume that this indexation rate is a certain 
proportion of the salary increase as this latter is highly positively correlated to the inflation. From 
the above definitions, uncertainty, represented by vector ((1 + wt ), affects boths assets and 

liabilities. As often in the literature (e.g. Kouwenberg m), we use a vector autoregressive model 
(VAR model) such that: 

ht = c + £lht—i + €t, €t ~ N (0, S), 

ht : = (In (1 + w t ) In (fi t ) ••• In (f fct ) ••• ln(f dt )) T , (2.2.6) 

t € 71 

where ht is a {(d+1) X 1} vector of continuously compounded rate, c the {(d+1) X 1} vector of 
coefficients, Q the {(d + 1) X (d + 1)} matrix of coefficients, et the {(d + 1) X 1} vector of error term 
and £ the {(d+1) x (d+1)} covariance matrix. The parameter estimation of this model requires 
time series analysis. For example in Kouwenberg [20], annual observations of the total asset returns 
and the general wage increase from 1956 to 1994 are used to estimate the coefficients of the VAR 
model. The resulting estimates will serve in constructing the scenario tree which constitutes the 
workhorse of multistage stochastic programs. 

2.3 The ALM problem 

The total cost of funding is the sum of regular (^ crt-iWt) and remedial (^ Zt) contributions over 
the studied period. In this study, we are looking for the investment stragtegy Ht, contribution 
rate crt and remedial contribution Zt for which the total expected cost of funding is minimized. 
The optimization is made under risk, legal, budget, regulatory and operating constraints. The 
constraints and objective of the ALM study will be presented in this section. 

We denote by symbol E* (x) the conditional expectation of random variable x with respect to 
the natural filtration J~t whereas P {E} denotes the occurence probability of event E. At each 
decision time t, the optimization problem consists in minimizing the total expected costs under the 
constraints considered in the following subsections. To simplify the notation, we omit the scenario 
index s. 

2.3.1 Risk constraints 

The pension fund wants to guarantee the participants a certain amount of pension. But the members 
also depend on the pension fund to actually provide for their needs in the future. Therefore, the 
safety of the portfolio is of paramount concern. This safety is translated in risk constraints. 

A pension fund has long term obligations, up to decades, and therefore, its planning horizon is large, 
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too. The main goal of an ALM is to find acceptable allocations which guarantee the solvency of the 
fund during the planning horizon. In general, solvency is measured by the funding ratio Iq (also 
called cover ratio) that we define for a given time t by 



Underfunding occurs when the funding ratio is less than one. The assertion Ft < 1 is equivalent to 
saying that the surplus at time t. i.e. A1 — Lt, is negative. When this occurs, the shortfall could be 
provided by the fund’s sponsor or any other external contribution. That is the remedial contribution 
as in Haneveld and al. |13j. Depending on how the random vector uit := ((1 + Wt) , £t) , t E 71, 
behaves, Ft may change over time. Therefore, the pension fund rebalances its assets portfolio and 
redefines its contribution rate in order to control the funding ratio. The higher T) is, the healthier 
the fund is. However, the decision maker would like to avoid as much as possible the changes in 
contribution rates. We will see in the model description that the parameters can be set in order to 
limit those variations. 

The long term objective of the pension fund consists in fullfilling both long and short (one year) 
term constraints. We define two types of funding ratio risk constraints in this paper. Their goal is to 
constrain the funding ratio to be larger, on average, than a predefined minimum 7,7 > 0. Namely, 
the expected shortfall E^_ 1 (At — yL/j) - , h > t, is required to be less than a certain amount (3t 
known at time t. Here, (a) - := max{—a,0} is the negative part of a E M. Also in order to simplify 
understanding, the expression expected shortfall is used to name E/,,_i (A* h — 7 Lh)~ ■ That is slightly 
different from its definition in actuarial science where 7 has to be equal to one. The one period risk 
constraint (OICC0) is expressed by 

Ei(A* +1 - 7 L m )- <fi t , teT 0 \{T}. (2.3.1) 

and for the multiperiod (MICC) approach, 

E fe _i (A* h - 7 L h )~ < fit, he Tt+x and t € 7o \ {T} (2.3.2) 


where fit and 7 are parameters defined by the pension fund. According to the short term approach 
at each decision time t, the expected shortfall over the following period should be smaller 
Notice that, at time t, the short term risk only controls the expected 
over the following one-year period. 


(2.3.1) 


than a certain amount fit. 
shortfall of (A * +1 — yLj+i) 

When we want to control the expected shortfall over the whole remaining period up to maturity, 


the risk constraint (2.3.2) is a good measure of long term risk (multiperiod). That is, at time t, 
equation (2.3.2) means that the one period expected shortfall E^-i (Ah — qL^) - should be smaller 


than fit at any future node with h E Tt+x- Equation (2.3.2) can be rewritten as 


max E fe _i (A* h - 7 L h ) < fi t , te%\ {T} 
h&Tt+i 


(2.3.3) 


meaning at time t, that, the highest one year expected shortfall, over the remaining periods to 
maturity, has to be smaller than the amount fit . Parameter fit is defined at time t by the pension 
fund and is a function of available informations at that time; e.g. fit := / (At, Lt ) = a At, 0 < a < 1 . 



shortfall, are named integrated chance constraints (ICC). They were proposed by Haneveld |17j 
as a quantitative alternative for chance constraints (CC). In section [3j both ICC and CC will be 
discussed more specifically. Successful applications of the ICC in ALM for pension fund can be 
found in Drijver ED and Haneveld and al. m- The authors assumed that fit := fi is unchanged 
over the studying period and in their numerical illustrations, the ICC is only applied to the first 


1 OICC (resp. MICC) stands for One period Integrated Chance Constraint (resp. Multiperiod Integrated Chance 


Constraint) which will be more clearly defined in section 3.2 
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stage. In such case, one can prove that the OICC and the MICC are equivalent. Instead, we remove 
that assumption in our work. Therefore, we define: 

0t ■■= aL t (2.3.4) 


where a , 0 < a < 1, denotes a scale free risk parameter. Until now, we are not aware of any 
implementation of the MICC (2.3.2) in a multistage framework. OICC is more relax than MICC. 
Obviously, the two constraints cannot be implemented in the same model at the same time. There, 
comes the other particularity of this paper: we analyze the multiperiod risk constraint and then 
measure how conservative it is comparing to the one period approach. 


2.3.2 Other constraints 


Risk constraints are important, but institutional and legal rules regarding pension fund operations in 
general are also relevant. As stated in Pflug and al. [28], institutional and legal rules are designed to 
restrict the risk of losses which would adversely affect pensionners. Hence, the following restrictions 
are integrated to the model. 

Firstly, the fund is not allowed to sell an asset that is not owned. This is the not short selling assets 
constraints and can be expressed by 

Hk,t > 0 , 

Bk,t > 0 , 

Sk,t >0 for k G /C and t £ To- 

The not short selling constraint goes with the not borrowing cash constraint expressed by 

C t >0, t E To- 


Secondly, at any time, the fund should dispose a minimum amount in cash in order to pay eventual 
claims such as death benefits or pensions. This can be called liquidity constraint and is formulated 
in our model as 

C t (1 + r f ) + E t {cr t W t+ i - Ben t+ 1 ) > 0 

which means that, on average, the cash allocation Ct at time t should be sufficient to cover the 
eventual negative value of the cash flow balance over period [ t, t + 1], Notice that the term liquidity 
constraint used here may have a different meaning in another context, e.g Fonseca and al. [12] in a 
macroeconomic framework. 

Thirdly, the fund is subject to portfolio constraint imposed by the legislator in order to keep a 
minimum control on its risk exposure. It consists on bounding the holding in asset k by setting 
upper and lower bounds, Uk and Ik respectively, on Hk,t- That is 


Ik-A-t C Hk : t C UkAt, k € 1C, t € To \ {T} . (2.3.5) 

For example in Switzerland the amount allocated to stocks should not exceed fifty percent of total 
wealth. In such case, l s tocks = 0 and u s t oc ks = 0.5 At and constraint (2.3.5) is equivalent to 

0 < H s t 0 cks,t C 0.5 At, t £ To \ {T} . 


These bounds are also applicable to cash Ct and we obtain 


IcAt <C t < u c A t , t £ To \ {T} . 


(2.3.6) 


Notice that, in equations (2.3.5) and (2.3.6), upper and lower bounds can also be time dependent. 
Finally in an asset allocation problem, dynamics and budget constraints, already defined in section 
|2.2| are inavoidable. If they were left out, the optimization program would be unbounded. The 
constraints presented in this subsection are common in any ALM stochastic programming imple- 
mention; see for e.g. Kusy and al. [21] . Carino and al. [7], Dempster and Consigli [9], Bogentoft 
and al. [6] and Dert |10| among others. 


OPP2 of April 18 th , 1984, Art 55-b, (As of January 1 st , 2012) 
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2.3.3 The optimization problem 

The ALM model is a dynamic decision making optimization tool to minimize the total expected cost 
under risk and operating constraints. Decisions are taken at the beginning of each one-year period. 
Accordingly, the ALM model is developped as a multiperiod decision problem, for which, we are 
asked to come up with an optimal asset allocation, contribution rate and remedial contribution at 
the beginning of each year. Moreover, penalty costs are assigned to the undesirable events: remedial 
contributions, and yearly absolute variation of contribution rates. All these components together 
constitute the objective function: 


min En 

H,cr,Z 


“T—1 T—2 

vt+1 (crtWt+i + \ z Z t +i) + ^t+iAA cr A crt Wt+i 

,t=0 t=0 


(2.3.7) 


where A crt := \ cry+i — crt | is the absolute variation of contribution rate from year t to t + 1, vt 
is the discount factor for a cash flow in year t. X z and Aa c are, respectively, penalty parameters for 
remedial contribution and absolute variation of contribution rate. The variables crt and A crt are 
bounded: 

cr l < crt < cr u and A cr < A crt < A cr , t € 7o \ {T} 

where cr l , A cr are the lower bounds and cr u , A cr the upper bounds of crt, A crt , respectively. The 
optimal decisions have to lead to a funding ratio greater than a certain minimum F (sometimes 
called target funding ratio) at the end of period of study T : 


At -f, 

Ft = — > F. 

Lt 

The entire ALM model, with objective and constraints, can be found in Appendix 1. An opti¬ 
mization program such as (2.3.7) is often referred to as a here and now problem. Uncertainty, 
characterized by ut = ((1 + Wt) , £t) , t E 71, is approached by scenarios. Therefore, we define Co 


with a finite number S of possible realizations u s := (d)f, ■ 
to t = T with relative probability p s 


, UJn 


€ S := {1,..., 5}, from t = 0 


The objective (|2.3.7|) is obviously linear as it can be rewritten as a linear combination of deci 
sion variables. 

\ and 


We can also notice that dynamics and constraints (except risk constraints which 
have to be rewritten in a linear form for the stochastic program) presented in sections 
are all linear in decision variables. If the risk constraints OICC (2.3.1) and MICC ( 2.3.2| ) were writ¬ 
ten in a linear form, the ALM problem would be a stochastic linear program (SLP), theoretically 
solvable by any SLP software depending on its size. In the next section, we will show how they can 
be turned into linear. The books Kail and Mayer m, Shapiro and al. |34| and Birge and Louveaux 
[5] provide good ressources to deal with such problems. When the size is big, resolution may require 
heuristic methods. Size is big means that number of asset classes is large or/and time horizon is 
long or/and number of scenarios is large. Decisions variables are Ht, Bt , St, Ct, crt and Zt for 
t E 7o; but only first stage values Hq, Bq, Sq, Cq, cro and Zq are crucial to the decision maker, 
since, almost surely, a true realization of the random data will be different from the set of generated 
scenarios. 

Stochastic programming (SP) is getting popular in ALMs. Its advantage lies in its ability to easily 
incorporate various types of constraints, Zenios and Bertola m- It has rooted with the work of 
Ziemba and al. m who showed, based on a 5-year period application to the Vancouver City Saving 
Credit Union, that SLP is theoretically and operationally superior to a corresponding deterministic 
linear programming (LP) model. The authors have proved that the effort required for the imple¬ 
mentation of ALM and its computational requirements are comparable to those of the deterministic 
model. Since then, SP in ALMs has been revisited by many other authors, see, for example, Ziemba 
and Mulvey m, Carino and al. [7J, Aro and Pennanen [3| and Zenios and Bertola m among 
others. 

By definition, the pension fund risk problem is often a shortfall problem. In such models, the 
relevant measure of risk for the firm is the expected amount (if any) by which goals are not met, 
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Carino and al. [fj. The model considered in this paper has a general DB ALM structure such as 
explored in Haneveld and al. [13] and Ziemba and al. m- Its main particularity consists in the 


integration of ICC by the way of OICC (2.3.1) and MICC (2.3.2). A successful implementation of 


constraints (2.3.1) in ALM for a DB fund can be found in Vlerk and al. m and Haneveld and al. 
[ 13] , In their works, the optimization problem is solved assuming that the parameter fit is constant: 


fit = fi. Furthermore, remedial contribution are provided only when funding ratio falls short in two 
consecutive years. Implementing this latter condition has lead to the use of binary variables. The 
authors proposed a heuristic solution to the problem. 

As we will explain in section 3.2, the parameter fi is not scale free. A certain value of fi does 
not have equivalent meaning for two different pension funds. It can be too low for a certain fund 
whereas too high for an other one. In addition, the pension fund actual situation should be taken 
into account. Our paper is an extension of Haneveld and al. m- As a novelty, we assume that the 
risk parameter fit vary with respect to time t and is defined as a proportion a of the actual level of 


liability at time t. see equality (2.3.4). Roughly speaking, on average, the total asset should cover 


a proportion of magnitude (1 — a) of liability at any time. In our model, remedial contribution can 
however be provided at anytime where solvency is in question, avoiding the use of binary variables, 
and indirectly, the need of heuristics. Penalty parameter punishes the abuse of remedial contri¬ 
butions. 

The main features of this study turns around the following points: 

• As in Haneveld and al. m where optimal decision is analyzed for different values of their risk 
parameter fi, we first measure the effect of our risk parameter a on the decisions Hq, ctq and 
Zq] this with respect to the OICC. In addition, for a fixed value a, the influence of the initial 
funding ratio is also explored. 


Secondly, as a safer alternative to the OICC, we propose the MICC (constraint (2.3.2)) and 


we then measure how hard it is, comparing to the OICC. In constraint (2.3.2), index h is a 


decision time index and we are not aware of any implementation of such constraint in ALM. 
The OICC considered in the first item is actually extended to a multiperiod risk constraint, 
reinforcing the long term aspect of the pension fund’s ALM. 

In the rest of this paper, OICC (resp. MICC) will stand for the one period (resp. multiperiod) ICC 
itself as well as the ALM model with the OICC (resp. MICC). 


3 Framework of the risk constraints 


The most important constraints, of course, deal with the goal of the pension fund: in all circum¬ 
stances keep a certain control on the funding ratio. This latter is expressed in terms of shortfall 
constraints which are of ICC type in this paper. Proposed by Haneveld H3, the ICC’s formulation 
direclty results from CC’s. That is why, in this chapter, we firstly introduce CC and how it leads 
to ICC. Secondly, ICC is discussed and we show how constraints (2.3.1) and (2.3.2) are related to 
it. We finish this section by proposing simple linear reformulations of (2.3.1) and (2.3.2). 

For the sake of clarity, we define the generic linear function G : xS-> 


such that 


G(X,u) := BX — D 

where X £ X is an d -vector of decision variables, X C is a polyhedral and closed set and 

oj := ( B,D ) : fi -» M m X R d X R m is a random parameter on the probability space (fi, T, P). 
The support of ui is defined as the smallest closed set S C R m X R rf X R m having the property 
P (w £ S) = 1. For i £ X := {1, • • • ,m}, the vector B is of dimension R m x such as B : = 
(Hi • • • Bi ■ ■ ■ H m ) T with Bi £ R rf whereas D := (Hi • • • D, • • • D, n ) T with D % £ R. As 
supposed in our SP model, we assume that uo = (B, D ) has a finite number S of possible realizations 
co s = ( B s , D s ), s £ S = {1,..., S} with respective probability p s . 
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3.1 Chance constraints 


Chance constraints (CC) models serve as tool for modeling risk and risk aversion in SPs. Let 0 be 
a m-dimensional vector of zeroes. Satisfying the constraint G (X, u) > 0 could lead to high costs or 
unfeasibilities. This equation refers to a finite system of m inequalities. Instead, if the distribution 
of u> is known, one can formulate the condition that the probability of G (X, uj) > 0 is sufficiently 
high, i.e. close enough to 1. That is 


P {G (X, w) > 0} > 1 — e 


(3.1.1) 


where the fixed parameter (1 — e) £ [0,1] is called probability level and is chosen by the decision 


maker in order to model the safety requirements. Equation (3.1.1) is the general form of chance 


(probabilistic) constraints and can be viewed as a compromise with the requirement of enforcing 
the constraint G (X, uj) > 0 for all values uj £ S of the uncertain data matrix. 


When m = 1, G(X, uj) := g(X,uj) is a scalar and equation (3.1.1) leads to 


P {g (X, u;) > 0} > 1 — e (3.1.2) 

with g : x E —> M. Equation ( |3.1 ~2j ) is known as individual CC. For m > 1, we obtain 

P {gi (X, uj) > 0, * G 1} > 1 - e, (3.1.3) 


called joint CC. 

Chance-constrained programs have been pionnered by Charnes and al. [8] in production planning. 
Since then, they have been extensively studied and have also been applied in many other areas 
such as telecommunication, finance, chemical processing and water ressources management. De¬ 
spite important theoretical progress and practical importance, there could be major problems with 
numerical processing of CCs, see Ahmed and Shapiro |2j and Nemirovski and Shapiro m- 
Especially when uj has a discrete distribution, Raike [30j introduces a mixed-integer reformulation 


of CC. Assuming m = 1, equation (3.1.2) is equivalent to 


^P S ■ l(g(x,u,‘)>o)(s) > 1 - e 


3=1 


where 1( 5 (x,uj s )> 0)( s ) = 1 if g{X,uj s ) > 0 and 0 otherwise. Now, we are able to write inequalities 
(3.1.2) in a mixed-integer program (MIP) formulation. We introduce binary variables 5 s , s E S. 


They play the role of indicator function: 5 s = 1 in scenario s if it holds that g(X,uj s ) < 0 and 
equals 0 otherwise. In terms of these additional decision variables, the CC can be written as linear 
inequalities 


g s ( X,uj s ) + 6 s M > 0, 
s 

Y^P s 5 s <e, 

S =1 

iGX, 5 s g {0,1} , 


s 6 5, 

(3.1.4) 

s £ S, 

(3.1.5) 

s £ S , 

(3.1.6) 


where M is a sufficiently large number. If 5 s = 0, then the constraint g (X, oP) > 0 corresponding 
to the realization s in the sample is enforced. On the other hand, if 5 s = 1, the constraint is satisfied 
for any candidate solution. The probability weighted average of these binary variables equals the 
risk of not meeting the condition g (X, u; s ) < 0 with the decision X, which should be at most e. 
This formulation is well known in SP and has first been applied to ALM for pension fund by Dert 
|10| . It also holds for the joint CC case where m > 1. In fact, {gj (X, u) >0, i G 1} is equivalent 
to 

min {gi (X, ca)} > 0 

i&l 
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and can also be writen as linear inequalities 


g!(X,u> a ) + 5!M> 0, 

i e /, 

s G S , 

(3.1.7) 

A s > Sf, 

o 

i <E /, 

S G S , 

(3.1.8) 

D 

c—1 


S G <S, 

(3.1.9) 

5— 1 

xeX, A s ,<5f e {0,1}, 

i e /, 

S G S . 

(3.1.10) 


Even with these linear settings (3.1.4) — (3.1.6) and (3.1.7) — (3.1.10), implementing this constraint 
with a reasonnable number of scenarios can be computationally challenging as the feasible set is 
obviously not linear, neither convex. That is due to the increase in complexity from MIP that arises 
from the introduction of at least one binary variable per each of the S scenarios. Efficient solution 
algorithms are proposed in chapter 4 of Kail and Mayer |15j , Luedtke [23], Luedtke and al. [24] . 
Tanner and Ntaimo [36], Prekopa and al. [29] and Ruszczynski [33] . 

Note that the CC, as described above, only considers the qualitative aspect of the risk, i.e. attention 
is only paid to whether the integrand is satisfied or not. A better approach can be to control the 
quantitative aspect of the fail, i.e the size of negative values of G s . That is often the case for pension 
funds where sponsors want to know approximatively how much they are willing to contribute in the 
following periods. Due to an idea of Haneveld binaries 6 s are dropped and the integrated 
chance constraint has been proposed. 


3.2 Integrated Chance Constraint 


The MIP constraints (3.1.7) to (3.1.10) are hardly implementable due to the integrality conditions 
in (3.1.10). For problems involving binary (or general integer) decision variables, a natural approach 
is to relax the integrality and solve the resulting relaxation, see Vlerk and al. m- If we relax the 
integrality constraints and substitute y s := 5 s M and /3 := aM , we obtain 


B s X + y s >D s , seS 

(3.2.1) 

5 

VI 

CO 

53b 

CO 

w 

(3.2.2) 

S=1 

y s > 0, s £ S 

(3.2.3) 

x ex, 

(3.2.4) 


where the parameter (5 is non-negative. By (3.2.1), for each s, the non-negative variable y s is not 
less than the shortfall ( B s x — D s )~, where ( a)~ := max{—a, 0} is the negative part of a £ M. The 


inequality (3.2.2) therefore puts an upper bound (3 on the expected shortfall. That is, the system 


(3.2.1) — (3.2.4) is equivalent to 


E (BX - D)~ = J2p S ( B ° X ~ DS )~ < P- 


(3.2.5) 


5=1 


Such constraint is called integrated chance constraint (ICC) and has been introduced by Haneveld 
m as an alternative to CC. However, Haneveld and al. [13], Vlerk and al. m and Drijver m 
have pionnered its application to ALM for pension funds and since then, it has been implemented 
in practice. 

By definition, the feasible set, defined by linear inequalities (3.2.1) — (3.2.4) is a polyedron (convex) 
as it contains only continuous decision variables, see Haneveld and Vlerk M- Thus, it can usually 
be solved efficiently using an appropriate software. Constraints (3.2.1) — (3.2.4) are very attractive 
from an algorithm point of view. Haneveld and Vlerk m propose a faster algorithm for big size 
problems. ICC is a good alternative to CC from different perspectives: 
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• Firstly, CC only measures the probability of shortage whereas ICC uses the probability distri¬ 
bution to measure the expected magnitude of the shortage. We can say that ICC takes into 
account both quantitative and qualitative aspects of the shortage whereas CC only considers 
its qualitative side. CC says only if there is underfunding or not and especially in practice, 
it could be important to limit the amount of remedial contributions the sponsor is willing to 
provide in years after. 


• Secondly, ICC and CC somehow ressemble, respectively, to the so-called conditional value- 
at-risk (CVaR) and value-at-risk (VaR). Conversely to CVaR which is known as coherent 
(Rockafellar and Uryasev [32]), it is well known that VaR is not a coherent risk measure as 
it does not fulfill the subadditivity condition. Therefore, ICC possesses more attractive risk 
properties than CC. To learn more about coherent risk measures, see Artzner and al. |4j. 

• We should also add that, if the risk aversion parameter is changed, the feasible region in case 
of ICC changes smoothly, while this region changes in a rough way in case of CC, Drijver HD- 

• Finally, we should admit that the parameter e of CCs is scale free, and corresponds to risk 
notion which is more familiar to pension fund managers. It is not the case for ICC. Our 
solution to this problem is to set (3 as a proportion a of liability. 


From now, and without loss of generality, we assume m > 1. 
rewritten as 


Therefore, equation (3.2.5) can be 


E {(BiX-Di)-, iel}</3 


which is the joint form of ICC, see Haneveld and Vlerk m- When index i is a decision stage index 
with conditional expectation at stage i, we obtain a multistage program and variable X becomes 
stage dependent (X t ). That is, at stage j £ / \ { m }: 


Ej (B i+ iXi - A+i) < /3j, i€ {j,j+ !,■•■ ,m-l} 


(3.2.6) 


which is equivalent to the MICC ( 2.3.2[ ) for / = 7o and Bh+iXh — -Dh+i = Ah + \ — ‘jLh+i. At time 
t, that is: 

Eh {A* h+1 - 7 Lh+i)" < ft, h G T t \ {T} . (3.2.7) 

The parameter j3t is then set at time t and will remain applicable until T. As decision is taken at 
each stage, the MICC inequality (2.3.2) shows a collection of inequality (3.2.7) going from t = 0 to 
t = T — 1. Similarly, when m = 1, one can proove that equation (3.2.5) leads to the OICC (2.3.1). 


3.3 OICC and MICC: Scenario tree interpretation 


Section 2.1 briefly explains our scenario tree model. We recall that the node (t, s ) corresponds to a 
certain scenario s at decision time t. To avoid anticipativity, we have to consider that many nodes 
(t. s ) might correspond graphically to the same thing on the scenario tree picture. For example in 
Figure 2.1.1] the nodes (1,1), (1,2), • • ■ , (1, 8 ) correspond graphically to the empty red cercle. At 
each node (f, s), the fund’s manager has to rebalance the asset portfolio and fix the contribution 
rate. These decisions are taken considering the actual scenario and possible future paths as well as 
the risk constraints. 


3.3.1 OICC 


In principle, considering a certain node ( t,s ), the OICC constraint (2.3.1) would be implemented 
as follows: 

(3-3.1) 


E m (A £! - 7-^t+i) := E Pi { A *U - 7-^t+i) 


<aL s t 


s£S 


where pf s stands for the conditional probability to reach node ^f+l,sj going from (t,s) and 

i , 

pf s = 0 for any scenario s of t + 1 not descending from (t, s ). As in Vlerk and al. [37], we include 
the linear inequality (3.3.1) in every subproblem (t,s ), t <T of our multistage recourse model. At 
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(t,s), they reflect the short-term risk constraint, stating that the expected funding shortfall over 
the following period (t + 1) is at most ctLf. In other words, on average, the pension fund should 
be able to cover the proportion (1 — a) of its total liability. The optimization problem will be more 
and more relax with the increase in a. 


3.3.2 MICC 

Considering the node (t, s ), the MICC constraint can be formulated in the following way: 

E h _ M (A* h s - 7 L£)- < aL s u h 6 Tt+i and t€%\ {T} (3.3.2) 

with 

E h . hs (. A* h s - 7 L s h )~ = £ P S h-i,s (V - 7 l()~ ■ 

s'es 


Under (3.3.2), at each node (t, s), decisions are taken such that the descending nodes’s one-period 
expected shorfall are smaller than aLf (defined at current node). Such constraint permits to have a 
certain control of the cover ratio over the whole remaining periods: [t + 1,T]; whereas (3.3.1) only 
covers one period: [t,t+ 1]. For example, at initial time t = 0, the minimum cost is defined such 
that the expected shortfall at any node in the tree (as descendant of the initial node) is smaller than 
/3 q = ocLq as in Haneveld and al. [13]: 


4 *S 

H+l 


< An t G % \ {T} , s e S. 


E p; 

s' SiS 

Futhermore, at each node (t, s) , t £ 71 \ {T} , s E S, we add the restriction: 


Epm 

s' £S 


A*S _ 

/i t+1 


< 


That is how we implement (3.3.2) at initial node. If we repeat the same procedure at each node of 
the tree, we can then propose a simpler SP reformulation: 


Proposition 3.1. Constraint (3.3.2) is equivalent to the following linear statement: 
At each node (t,s ), t < T, s £ S 


s' G.S 


Pu 


4* s 
^t +1 


< 


min aLf;. 

o <t'<t 


(3.3.3) 


That is, at a given node (t,s) ,t < T,s € S , the expected shortfall over the next period should 
be less or equal to the smallest value of aLf calculated over the preceding nodes (t 1 , s ), t' < t. This 
is based on the fact that, in the multiperiod framework, decision taken at node (t, s) is influenced by 


the history of cuf up to time t. in particular /3* at preceding nodes. Inequality (3.3.3) is linear and 
describes a polyhedral set. The proof of proposition |3 .1 1 is straightforward when we go backward in 
time starting from nodes (T — 1, s), see Appendix 2 for an example based sketch of proof. At each 
node (t,s), as we know the history of j3f up to time t, one can determine the smallest /3|,, t' < t. 
Therefore, implementation of MICC consists in including the linear constraint (3.3.3) at each node 
(t,s). 


4 Numerical illustrations 

This section contains computational results for the SP model. Let’s recall that we are dealing with 
a DB pension fund whose objective is to minimize the total expected costs under constraints. The 
study will focus on risk constraints which are of ICC type. Firstly, based on the OICC, the effects 
of risk parameter and cover ratio on the optimal decisions are analyzed. Prima facie, the MICC 
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Table 1: Data on the asset classes 


Asset classes 

k 

Ik 

U'k 

c 

Initial Investments 

Cash 

— 

0 

l 

0 

4’950 

Deposits 

1 

0 

0.5 

0.00150 

16’500 

Bonds 

2 

0.1 

1 

0.00150 

38’500 

Real estate 

3 

0 

0.30 

0.00425 

17’600 

Stocks 

4 

0 

0.50 

0.00425 

32’450 


Table 2: Values of the other deterministic parameters 


\ z = 350 A Acr = l r f = 0.008 v t = (1 +r f )~ t 
A cr = —0.08 cr l = —0.08 A cr = 0.05 cr u = 0.3 
17= llO'OOO 7 = 1.05 F = 1.05 


appears to be a safer and more restrictive than OICC. Based on the same analysis as before, the 
cost of conservativeness is subsequently measured. 

For this study, consider a hypothetical large pension fund which may invest into d = 4 classes of 
asset ordered by level of risk: 

1. Deposits, 

2. Bonds, 

3. Real estate, 

4. Stocks. 

We are aware of the fact that the number of assets is often much larger in practice. That said, 
only four classes of assets are considered here in order to reduce the complexity of the model. After 
investing in these asset classes, the rest is held in cash. The deterministic properties of asset classes 
are described in Table [lj Investment limits are defined with respect to practical rules of liquidity 
and diversification; transaction costs are taken from Haneveld [13] with c s = c B = c, whereas the 
initial investments are defined considering general statistics of pension fund’s assets allocation in 
Switzerland, see Towers Watson (38] (where we assume that "real estate" corresponds to "other 
assets"). The portfolio constraints are defined in term of proportion and all amounts are assumed 
to be in thousands of Swiss francs. The values of the other deterministic parameters are shown in 
Table [2] 

The time horizon T = 5 is split into five periods of one year each. Consequently, the considered 
ALM model has five stages, allowing for decisions at t = 0 (now) up to time t = 4. The random 
vector oot follows a VAR process, approximated in our case by a multistage scenario tree. In the 
following considerations, we firstly present the descriptive statistics of our model and the numerical 
results, obtained from our study, are discussed on a second hand. 

4.1 Scenarios 

The implementation of the scenario tree requires a careful specification of the VAR process. For this 
purpose, we use the estimation results obtained in Kouwenberg [20]. More specifically, the author 
estimates this process based on annual observations of the total asset returns and the general wage 
increase from 1956 to 1994. Table [3] displays descriptive statistics of the time series whereas Table 
[4] shows the estimated correlation matrix of the residuals. Future returns for financial planning 
models can be constructed by sampling from the error distribution of the VAR model and applying 
the estimated equations of Table [5] We refer to Kouwenberg [20] for further details on this model 
estimation and for building the tree as well. For this purpose, we specify a branching structure of 
1 — 10 — 6 — 6 — 4 — 4. This scenario tree has one initial node at time 0 and 10 succeeding nodes 
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Table 3: Statistics, time series 1956-1994, Kouwenberg 


Assets 

Statistics 

Mean 

S.D. 

Skewness 

Kurtosis 

Wages 

0.061 

0.044 

0.434 

2.169 

Deposits 

0.055 

0.025 

0.286 

2.430 

Bonds 

0.061 

0.063 

0.247 

3.131 

Real estate 

0.081 

0.112 

-0.492 

7.027 

Stocks 

0.102 

0.170 

0.096 

2.492 


Table 4: Residual correlations of VAR-model, Kouwenberg [20] 


Assets 

Wages 

Deposits 

Bonds 

Real estate 

Stocks 

Wages 

1 





Deposits 

0.227 

1 




Bonds 

-0.152 

-0.268 

1 



Real estate 

-0.008 

-0.179 

0.343 

1 


Stocks 

-0.389 

-0.516 

0.383 

0.331 

1 


at time 1, ■ ■ ■, resulting inl0x6x6x4x4 = 5 = 5760 path from 0 to 5, each with probability 
I Y = 5760' 

4.2 Numerical results 

This section presents the outputs of our study. All numerical results were implemented using the 
solver CPLEX in the mathematical programming language AMPL. The ALM models are formulated 
as large LP-problems with 616 7 321 variables. In the model with the OICC, there are 995 7 347 
constraints and 3 / 041 / 032 nonzeros in the constraint matrix whereas they are respectively l / 002 / 317 
and 3 / 105 , 602 in the MICC. On average, the solution times are 381 seconds and 448 seconds, 
respectively, for OICC and MICC. 

As a result of the ALM analysis, we are supposed to provide the first stage optimal decisions: a 
contribution rate, a remedial contribution and asset allocation that minimize the total cost. In the 
first part of this section, we analyze the effects of the risk parameter a and the initial funding ratio 
Fq on the optimal decision. The optimization is made with respect to OICC. The values of a ranges 
from 0 to 0.085 whereas the initial funding ratio Fq = vary from 0.5 to 1.5. In order to vary Fq, 
we change the initial liability Lq accordingly, as Aq is specified from Table [ 2 ] In the second part, 
we compare the OICC to the MICC. 

4.2.1 OICC 


In order to analyze the impact of the risk parameter a, we fix the value of the initial cover ra¬ 
tio. That is: 

An 

Lo := 120 7 000 ^ Fq = — = 0.9166, under covered . 

-00 

In what follows, the letter O at the top of a symbol stands for OICC whereas M is related to MICC. 
Figure 4.2.1 shows the evolution of the contribution rate cro- The value of cro is particularly high 
as the institution is underfunded. We osberve that when a < a® := 0.025 (O at the top stands for 
OICC.), the contribution rate is at its maximum: cr u = 0.3 as specified in Table [ 2 ] From 0.025, 
cro decreases linearly until it reaches the value 0.27 at a = a° := 0.04, and remains unchanged 
thereafter. According to the objective (2.3.7), the total cost is composed of the total contribution 
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Table 5: Coefficient of the VAR model, Kouwenberg [20] 


In (1 + wages*) 

= 0.018 

+ 

0.693 In (l + wages* __*) + e** 

erg* = 0.030 


(2.058) 


(5.789) 


In (1 + deposits*) 

= 0.020 

+ 

0.644 In (l + deposits*.*) + * 

o’ 2 * = 0.017 


(2.865) 


(5.448) 


In (1 + bonds*) 

= 0.058 

+ 

e 3t 

°3,* = 0.060 


(6.241) 




In (1 + real estate*) 

= 0.072 

+ 

1 

05 ,* = 0.112 


(4.146) 




In (1 + stocks*) 

= 0.086 

+ 

^6 1 

<Js,t = 0.159 


(3.454) 






Figure 4.2.1: OICC: Contribution rate at t — 0 Figure 4.2.2: OICC: Initial cost allocation in 

as function of a function of a 


and of the total remedial contribution, these, over the period under study. Remedial contribution 
should be seen as an external financial support which may come from the sponsor of the pension 
fund. Figure |4 .2.2 1 displays the allocation of the total costs into the two types of contribution. The 
proportion of remedial contribution linealy decreases from 9% at a = 0 to reach 0% at a° = 0.025 
and stays constant for a > a®. Indeed, Figures 4.2.1 and 4.2.2 help understand how the ALM 


model parameters have been defined. It is conventional to assume that, from a certain level of risk 
and for a fix cover ratio, the sponsor will no more provide any financial support to the fund. In this 
model, the penalty parameter has been set such that the total remedial contribution is zero for 
a > a°. Consequently, cro decreases from a = a]p. It remains equal to 0.27 for a > a° due to the 
target cover ratio at maturity: F§ > F = 1.05. 

Figure 4.2.3 describes the optimal asset allocation for different values of a. Assets are ordered with 
respect to their level of risk. For small values of the risk parameter, the proportion of riskier assets 
(stocks and real estate) tends to increase with an increase of a. When it approaches a as the 
remedial contribution is already low, the decision maker starts reducing the risk exposition of its 
portfolio. However, the proportion of bonds is increased in order to improve the performance of 
the asset portfolio. The risk exposition is then progressively reduced until a = or, from which, it 
remains unchanged thereafter. The value a° can be seen as the smallest value of a, above what, 
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Figure 4.2.3: OICC: Asset allocation at t = 0 in function of a 


the OICC influence is no more significant, i.e. contribution rate, remedial contribution and asset 
allocation stay constant. 



Figure 4.2.4: OICC: Contribution rate at t — 0 Figure 4.2.5: OICC: Initial cost allocation in 

in function of F 0 function of F 0 


Next we consider the effect of the initial funding ratio on the first stage optimal decision. We 
vary the value of Lq so that the initial funding ratio Fq lies between 0.5 and 1.5, and we assume 
a = 0.035. Figure 4.2.4 displays the evolution of the contribution rate cro whereas Figure 4.2.5 shows 
how the total cost is distributed into regular and remedial contributions. As explained earlier, it is 
conventional to assume that, above a certain funding ratio Fq* , the remedial contribution is zero 
and contribution rate decreases as well. According to Figures 4.2.4 and 4.2.5 the ALM model is set 
such that Fq* := 0.9. 

Figure 4.2.6 describes the behaviour of the first stage optimal asset allocation with respect to Fq. 
When Fq < Fq , the optimal asset allocation is stable: approximatively 30% in riskier assets. From 
F 0 = Fq* , the cover ratio is high enough to no more obtain remedial contribution and to reduce the 
contribution rate. However, the decision maker has to act in a riskier way in order to meet pension 
fund liabilities. As a result, the risk exposition increases up to 50% at F$ = 1.275. An important 
target of our model is to guarantee a funding ratio F$ > F by minimizing the total cost and risk 
level. For larger values of Fq, with a higher chance to fulfill the condition F$ > F, the contribution 
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Figure 4.2.6: OICC: Asset allocation at t — 0 in function of F 0 



rate and the risk exposition (i.e. proportion of higher risk assets) decrease. We recall that the 
objective is not to maximize the wealth, but to minimize the total cost. Thus, the wealthier the 
fund is, the more prudent the allocation will be. 


4.2.2 MICC 


In this section, we firstly present the results of the ALM model with the MICC, and secondly, we 
compare with the OICC. For the MICC analysis, assumptions are similar to the ones made for the 
model with OICC. Figures 4.2.7 to |4.2.12 display the results of the analysis. The effect of the risk 
parameter a is measured in Figures |4.2.7[ |4~2.8| and [4.2.9 whereas Figures [4.2.10 4.2.11 and 4.2.12 


analyze the initial funding ratio impact. Although the first stage optimal decisions are different, 


they behave similarly. The parameters , 
slightly differ: 


of := 0.027; 


a M and Fq 1 * 


a M ■= 0.07; 


(respectively a|p, a° and Fq* for OICC) 


F 0 m * := 0.9. 



■ Contribution □ Remedial 





alpha 



Figure 4.2.7: MICC: Contribution rate at t — 0 Figure 4.2.8: MICC: Initial cost allocation in 

in function of a function of a 


The MICC approach is basically more rigid than OICC. That is why the above parameters are 
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Figure 4.2.9: MICC: Asset allocation at t — 0 in function of a 
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Figure 4.2.10: MICC: Contribution rate at t = 0 Figure 4.2.11: MICC: Initial cost allocation in 

in function of F 0 function of F 0 


greater or equal to their analogues in OICC. Notice that, according to Figures 4.2.3| and |4.2.9| and 
for a > oP in OICC and a > a AI in MICC, the asset allocations are exactly the same; showing that 
when a is above a° (resp. a M ), the OICC (resp. MICC) has no influence on the initial ALM model. 
In general, the optimal decisions obtained from OICC and MICC slightly differ. For example, when 
a = 0.05 and Fq = 0.9166, the first stage optimal decision of the MICC is: 


H 0 := (0.41 0.59 0 0) 


erg = 0.279 


whereas for the OICC: 


Hq := (0.26 0.66 0.08 0) T ; cr 0 = 0.270; 


and 


and 


Zn — 0 


Zn — 0. 


Asset allocations are calculated as percentage of the total asset. According to the above example, 
the decisions related to the OICC approach are riskier than the ones of MICC, especially regarding 
the asset allocation. In what follows, we will try to quantify the cost of this risk reduction. For 
a pension fund, this can be done by measuring the difference in term of the total cost (regular 


contribution + remedial contribution). Hence, Figure 4.2.13 compares the total costs of OICC and 
MICC whereas Figure [4.2.14| displays the contribution rate difference, all this with respect to a. 
When a < a !p = 0.025 or a > a R1 = 0.07, the contribution rate and total cost are equal for both 
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Figure 4.2.12: MICC: Asset allocation in function of F 0 
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Figure 4.2.13: Comparison of OICC and MICC 
in function of a: Total cost 
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Figure 4.2.15: Comparison of OICC and MICC 
in function of F 0 : Total cost 
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Figure 4.2.14: Comparison of OICC and MICC 
in function of a: Contribution rate cro 
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Figure 4.2.16: Comparison of OICC and MICC 
in function of F 0 : Contribution rate cro 
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models. For a jp < a < a M , OICC and MICC slightly differ, i.e. MICC costs more for a maximum 
variation of 2 , 000 (less than 2% of total asset) and 1.5%, respectively, for the total cost and the 
contribution rate. Consequently, although being more conservative, the MICC is preferable the 
OICC for the two following reasons: 


• it is safer, and 

• the cost of this safety is not significant: less than 2% of total asset. 

Moreover, Figures [4. 2.15 and 4.2.16 which compare the effect of Fq on the OICC and on the MICC, 
confirm that result. 


5 Conclusion 

In this paper, we considered the effects of integrated chance constraints (ICC) on an ALM model 
for a defined benefit pension fund. ICC are appropriate for modeling risk constraints, in particular 
when a quantitative risk measure is preferable. At each decision time, they put an upper bound on 
the one period expected shortfalls. The upper limit considered is a fixed proportion of total liability 
making the risk parameter scale free as well as time dependent. ICC are very attractive from a com¬ 
putational perspective. Their impact on the ALM model is analyzed through the implementation 
of a multistage stochastic linear program. We defined two types of ICC: the one period integrated 
chance constraint (OICC) and the multistage integrated chance constraint (MICC). 

The first step of the work consisted on describing the influences of the ICCs on the time 0 opti¬ 
mal decisions. This impact is measured by the way of the risk parameter and of the initial cover 
ratio. As one could expect, the risk exposition of the optimal portfolio tends to increase with the 
risk parameter until this latter reaches a certain value, from which, the allocation remains stable 
thereafter. On the other hand, for reasonable values of the initial cover ratio, the risk exposition 
is increasing. Above a certain value, as there is a lower chance of not meeting the target funding 
ratio condition and in order to guarantee the benefit payments at a low cost, the risk exposition 
decreases. 

MICC is basically more restrictive and more cautious than OICC. Secondly, the cost generated by 
the rise in security with the MICC is quantified. Although the optimal decisions from the OICC 
and the MICC are not the same, the total costs are very close, showing that the MICC is definitely 
a better approach. 

However, the obtained results are subject to discussion since they mainly rely on the scenario tree. 
In further considerations, it would be interesting to analyse the result stability regarding the sce¬ 
nario tree. 
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APPENDICES 

Appendix 1: The ALM program description 

We start this section by defining indices, variables and parameters of the model. Secondly, the ALM 

model with the objective and the constraints are also displayed. 
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Indices 


t 


time index, t = 0,1, • • • , T 
scenarios index, s = 1 , • ■ ■ , S 
index of asset classes, k = 1 , • • • ,d 


s 
k 

Decision variables 

Zf remedial contribution by the sponsor at time t in scenario s 

Cf total cash amount at the beginning of year t in scenario s 

Hf t value of the investments hold in asset class k, at the beginning of year t and in scenario s 
Bf t value of the asset class k , bought at the beginning of year t and in scenario s 

SI t value of the asset class k , sold at the beginning of year t and in scenario s 

erf contribution rate for year t + 1 in scenario s 

Af total asset value at time t in scenario s 

Af s total asset value just before the asset allocation and the remedial contribution at time t in scenario s 

A s variation (increase or decrease) of contribution rate from year t to t + 1 in scenario s 

Random parameters 

rf. t random rate of return on asset class k over year t in scenario s and t := 1 + r^t 

W t ? random total wages of active participants in year t in scenario s 

Lf random value of liabilities at time t in scenario s 

Benf random total benefit payments to active participants in year t in scenario s 

Deterministic parameters 

T time horizon 

S number of scenarios 

d number of asset classes 

a risk parameter defined by either the sponsor or the regulator 

cf proportional transaction cost for purchasing an asset class k 

o 

c% proportional transaction cost for purchasing an asset class k 
Ik lower bound on the proportion of asset class kin the total asset portfolio 
Uk upper bound on the proportion of asset class kin the total asset portfolio 
l c lower bound on the proportion of cash Adn the total asset portfolio 
u c upper bound on the proportion of cash Adn the total asset portfolio 
cr l lower bound on the contribution rate 
cr u upper bound on the contribution rate 

A cr lower bound on the yearly absolute variation of the contribution rate 
A cr upper bound on the yearly absolute variation of the contribution rate 
F target funding ratio 

7 lower bound on the funding ratio 

rf risk free interest rate 

penalty parameter for remedial contribution 
A Acr penalty parameter for absolute variation of contribution rate 
Hk,o value of the initial allocation in asset class k 

Co initial cash amount 
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Objective 

The objective of the model is to determine the asset allocation, contribution rate and remedial 
contributions that minimize the total cost defined as follows 

T-l T -2 

min E 0 Yf v t +i (cr t W t +1 + X z Z t+ i ) + Yf ut + iAA cr A crt W t+ i 

H,cr,Z — 4 *■—' 

U=0 t =0 J 

under the following constraints. In the model description, anytime we use indices s and/or k is 
equivalent to saying for any s = 1, • • • , S and/or for any k = 1, • • • , d. 

Budget constraints and total value of the assets 

d d d 

U = H kfi + Co + Zq — yy (cf + c k S k ,o) = Bfco + Co 

/c=l fc=l fc=l 

d 

A s t = yy + Cf_, (1 + rf) + cr^_ 1 Wf — Ben^ = Aj? 

k= 1 

d d 

At = Y + CU (1 + rf) + crUWf - Ben? + Zf - Y (cfBf t + cfSf t ) ; t = 1, • • ■ , T - 1 

fc=l fc=l 

d d 

At = at+ zt- Y^k B u+4 s u) = Y H it +c t' 

fc=i 

Asset classes dynamics 


k= l 


74,0 — + 4c,o — *‘4,0 

frs _ fS rrs . rjS cs + _ 1 T 1 

_ £k,t n k,t-l n k,t ~ ‘-’fed ) 1 ~ T ' ' ' > 1 


Cash dynamics 


Co — Co + Zq — yy (i+ cf) -B/^o+ y~f (1 _ cf) s k) $ 

k= 1 /c=l 

d d 

C s t = CU (1 + r f ) + crUWf - Ben? + Z\ - £ (l + cf) Bf t + £ (l - cf) ; 

fc=l /c=l 

IVof short selling assets and not borrowing cash constraints 

Htt> 0; Sf t >0; £4 > 0 ; Cf > 0 ; t = 0,---,T 


t = !,••• , T 


Liquidity constraints 


Cf (1 + r f ) + E t)S (cr?W t+ i - Ben t+ 1 ) >0; t = 0, • • • , T - 1 


Portfolio constraints 


l k A s t < H kt < u k A s t ; f = 0, ■ ■ ■ , T — 1 
l c A s t < Cf < u c A s t ; t = 0, • • • , T - 1 


Constraints on contribution rates 

cr l < cr s t < cr u and A cr < A? rf < A cr ; f = 0, • • • , T — 1 

The decision variables are subject to the non-anticipativity constraints. The integrated chance 
constraints defined in section [2. 3. 1| control the risk-level of the model. They also have to be included 


in the model as explained in section 3.3 
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t — 0 event t — 1 event t — 2 event t = 3 
Figure 5.0.17: A scenario tree with 40 scenarios and 66 nodes. 


Appendix 2: Proposition 3.1 > an example based sketch of proof 


We consider the event tree of Figure |5.0.17 with a time horizon T = 3 and a branching structure 
of 1 — 5 — 4 — 2, leading to5 = 5x4x2 = 40 scenarios. A node is a possible outcome of the 
stochastic event at a given time. The starting and ending nodes of the tree are round whereas the 
others are rectangular. At each t £ {2, 3}, the rectangular nodes are name according to time and 
following a top-down alphabetic order. For example, the rectangle (l,a) describe the outcome at 
the first node of time 1. At each node, the economical values such as total asset, total liability and 
expected shortfall can be determined. For sake of clarity, we recall that, due to non-anticipativity, 
the node (1, a) is equivalent to the node (1, s '), s' £ {1, • • • , 8} as described before, the node (1, b) 
is equivalent to the node (1, s '), s' £ {9, • ■ ■ , 16} and so on for the other nodes. The other notations 
used here are similar to the ones in the paper. We define 


A (M) := p U ( A 


1*S 

h+i 


s'eS 


where pf s stands for the conditional probability to reach node ^t+l,s^ going from (t,s) and 


3.3.2 


pf s = 0 for any scenario s' of t + 1 not descending from ( t , s). The MICC defined in equation 
is then 

{Afe-gs, h £ 77-i-i } £ Pt,si t £ { 0 , 1, 2 } . 

The value A,s ■= oLf is the ICC upper limit computed at time t in scenario s. According to this 
set of constraints, 


from initial node at t := 0 


h £ (1, 2, 3} A 0 , s < An M,s < An and A 2 , s < A), s ^ $■ 


(5.0.1) 


At t := 1, 


from node (1, a), h £ {2, 3} =n Ai jS < Al,s, and A 2)S < Al.sj s £ {1, 2, • • • , 8} 
from node (1,6), /i £ {2, 3} =n Ai jS < /3i,s> and A 2]S < Al.sj s £ (9, 10, • • • , 16} 


from node (1, e), h £ {2, 3} =n Ai jS < Al,s> and A 2jS < Al,s> s £ {33, 34, ■ 


40} 
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At t := 2 


from node (2 ,a), h = 3 =>■ and A 2 , s < 02,s, s G {1, 2} 
from node (2, 6), h = 3 => and A 2 )S < /? 2 ,s, s G {3, 4} 


from node (2,t), /t = 3 =>• and A 2 jS < /? 2 ,s, s G {39, 40} 
Therefore, we obtain for any s G 5 that 


at time 0, Ao, s < A) 

at time 1, A l; . s < A), A ljS < /3i, s 

at time 2, A 2 , s < A). A 2jS < /3i, s , A 2 , s < P\, s 


^-0,s A A) 

Ai jS < min {An /3i, s > 
^ A 2)S < min {/3q, Al )S , /? 2 ,s} • 


The obtained result leads obviously to the proposition 3.1 in the paper. Considering such an 
example is therefore without loss of generality. For an other event tree with different time horizon 
and branching structure, the proposition can be prooved using the same procedure. 
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